OCHA-PP-273, KEK-TH-1149, VPI-IPNAS-07-06 

Constraints on New Physics from 
Long Baseline Neutrino Oscillation Experiments 

Minako Honda/'Q Yee Kao,^'Q Naotoshi Okamuraj^'Iil 
Alexey Pronin,^'!^ and Tatsu Takeuchi^'Q 

^Physics Department, Ochanomizu University, Tokyo 112-8610, Japan 

^ ! ^Institute for Particle, Nuclear, and Astronomical Sciences, 

O 



o 



Physics Department, Virginia Tech, Blacksburg VA 24061, USA 



W). ^KEK Theory Division, Tsukuba 305-0801, Japan 

^3 ■ 

<^ : (Dated: July 31, 2007) 

Abstract 



New physics beyond the Standard Model can lead to extra matter effects on neutrino oscillation if 
Q^' the new interactions distinguish among the three flavors of neutrino. In a previous paper we 

argued that a long-baseline neutrino oscillation experiment in which the Fermilab-NUMI beam in 



31 is aimed at the planned Hyper-Kamiokande detector 



a 



would be capable 



■ its high-energy mode 

> : 

lO of constraining the size of those extra effects, provided the vacuum value of sin 2^23 is not too close 

' to one. In this paper, we discuss how such a constraint would translate into limits on the coupling 

■ constants and masses of new particles in various models. The models we consider are: models with 

o . 

, generation distinguishing Z s such as topcolor assisted technicolor, models containing various types 



of leptoquarks, R-parity violating SUSY, and extended Higgs sector models. In several cases, we 



^ ■ find that the limits thus obtained could be competitive with those expected from direct searches 

H : 

at the LHC. In the event that any of the particles discussed here are discovered at the LHC, then 
the observation, or non-observation, of their matter effects could help in identifying what type of 
particle had been observed. 
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I. INTRODUCTION 



When considering matter effects on neutrino oscillation, it is customary to consider only 
the ly-exchange interaction of the Vf. with the electrons in matter. However, if new inter- 
actions beyond the Standard Model (SM) that distinguish among the three generations of 
neutrinos exist, they can lead to extra matter effects via radiative corrections to the Zvu 
vertex, which effectively violate neutral current universality, or via the direct exchange of 
new particles between the neutrinos and matter particles 

Many models of physics beyond the SM introduce interactions which distinguish among 



9l models 



generations: gauged La — Lp [5[ and gauged B — aL^ — I3L^ — 'jL^ 
introduce Z's and Higgs sectors which distinguish among the three generations of leptons; 
topcolor assisted technicolor treats the third generation differently from the first two to 

i ngs in supersymmetric models 



explain the large top mass 
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12| : R-parity violatin g co up 



couple fermions/sfermions from different generations 
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15| 



The effective Hamiltonian that governs neutrino oscillation in the presence of neutral- 



current lepton universal 
differently, is given by 



ity violation, or new physics that couples to the different generations 
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In this expression, U is the MNS matrix 16| . 

a = 2EVcc , Vcc = v^G^iV^ = iV, 



4M^ ' 



(2) 



is the usual matter effect due to PF-exchange between i/g and the electrons [17], and b^,, b^, 
br are the extra matter effects which we assume to be flavor diagonal and non-equal. The 
matter effect terms in this Hamiltonian can always be written as 
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The unit matrix term does not contribute to neutrino oscillation so it can be dropped. We 
define the parameter ^ as 

hr - _ ^ 



(4) 



Then, the effective Hamiltonian can be written as 
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where we have absorbed the extra 6-terms in the (1, 1) element into a. 

The extra ^-dependent contribution in Eq. can manifest itself when a > \S'm'^^\ {i.e. 
E >10 GeV for typical matter densities in the Earth) in the z/^ and z/^ survival probabilities 



as 



ll 



P{ly^. ^ z/^) ^ 1 - sin^ 2^23 - 
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(7) 



and the CP violating phase 6 has been set to zero. As is evident from these expressions, the 
small shift due to C, will be invisible if the value of sin^ 26'23 is too close to one. However, if 
the value of sin^ 26*23 is as low as sin^ 26*23 = 0.92 (the current 90% lower bound [18|), and 
if ^ is as large as ^ = 0.025 (the central value from CHARM/CHARM II (igj), then the 



shift in the survival probability at the first oscillation dip can be as large as ~ 40%. If the 

™™, a.. 

toward the planned Hyper-Kamiokande detector [3[ in Kamioka, Japan (baseline 9120 km), 
such a shift would be visible after just one year of data taking, assuming a Mega-ton fiducial 
volume and 100% efficiency. The absence of any shift after 5 years of data taking would 
constrain ^ to [l| 

lel < ^0 = 0.005 , (8) 



at the 99% confidence level. 

In this paper, we look at how this potential limit on C, would translate into constraints 
on new physics, in particular, on the couplings and masses of new particles. As mentioned 
above, the models must be those that distinguish among different generations. We consider 
the following four classes of models: 

1. Models with a generation distinguishing Z' boson. This class includes gauged L^ — L^, 
gauged Le — L^, gauged B — aL^ — PL^ — •jL^-, and topcolor assisted technicolor. 

2. Models with leptoquarks (scalar and vector). This class includes various Grand Uni- 
fication Theory (GUT) models and extended technicolor (ETC). 

3. The Supersymmetric Standard Model with R-parity violation. 

4. Extended Higgs models. This class includes the Babu model, the Zee model, and 
various models with triplet Higgs, as well as the generation distinguishing Z' models 
listed above. 

These classes will be discussed one by one in sections II through V. The constraints on these 
models will be compared with existing ones from LEP/SLD, the Tevatron, and other low 
energy experiments, and with those expected from direct searches for the new particles at 
the LHC. Concluding remarks will be presented in section VI. 

II. MODELS WITH AN EXTRA Z' BOSON 

Z' generically refers to any electrically neutral gauge boson corresponding to a flavor- 
diagonal generator of some new gauge group. Here, we are interested in models in which the 
Z' couples differently to different generations. The models we will consider are (A) gauged 
Le — L^ and Le — Lr, (B) gauged B — aLe — PL^ — jLr, with a + P + 'j = 3, and (C) topcolor 
assisted technicolor. 

A. Gauged Lg — and Lg — L^- 

In Ref. 0, it was pointed out that the charges Le — L^, L^ — Lt-, and L^ — Lt are anomaly 
free within the particle content of the Standard Model, and therefore can be gauged. Models 
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FIG. 1: Diagrams that contribute to neutrino oscillation matter effects in (a) the gauged Lg — 
model, and (b) the gauged Lg — Lj. model. 



with these symmetries are recently receiving renewed attention in attempts to explain the 
large mixing angles observed in the neutrino sector Of these, gauged Le—L^ and Lg—Lr 
affect neutrino oscillation in matter. These models necessarily possess a Higgs sector which 
also distinguishes among different lepton generations 2l| , but we will only consider the effect 
of the the extra gauge boson in this section and relegate the effect of the Higgs sector to a 
more generic discussion in section V. 

The interaction Lagrangian for gauged — Li {i = fi or r) is given by 



c = 9z' [n^e - h^^ + ^i^^eL - vEYi^iL ) ^; . 



(9) 



The diagrams that affect neutrino propagation in matter are shown in Fig. [H (The exchange 
of the Z' between the Ue and the electrons do not lead to new matter effects.) The forward 
scattering amphtude of the left-handed neutrino uil {i = /i, t) is 



iM 



{Wz'){-Wz') ^ehY^eL Wei) I ttt ) (^1 ^7''e |e) 



(10) 



The electrons in matter are non-relativistic, so only the time-like components of the currents 
need to be considered. Replacing (e|e7°e|e) = {e\e^e\e) with Ne, the number density of 
electrons in matter, and (z/£l| z7^7°z/fL lun) = \ ^Il^£l Wil) with (f)l^4>ue, where is the 
wave function of the left-handed neutrino u^, we obtain 
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(11) 



Therefore, the effective potential felt by the neutrinos as they traverse matter can be iden- 
tified as 



V,, 
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UA 


7.6 


[22] 


OPAL 


10.6 


12.7 


8.6 


[23] 


DELPHI 


13.9 


12.2 


15.8 


m 


ALEPH 


12.5 


10.5 


12.8 


[2_5] 



TABLE I: The 95% confidence level lower bounds on the compositeness scale A^ (TeV) from 
leptonic LEP/LEP2 data. Dividing by -v/ivr converts these limits to those on {Mz' / Qz')- 



The effective ^'s for the Lg — and Lg — Lr cases are 



Vcc 



+ 



+4 



-4 



1 



1 



V2Gf \Mz' 



9z' 
Mz' 
9z' 



(13) 



Vcc ' {gVMl 

Ignoring potential contributions from the Higgs sector, a bound on ^ of |^| < = 0.005 
from Eq. ([8]) translates into: 

Mz' 



> 



1 



3500 GeV 



(14) 



9z' V V2Gf^o 
for both the L^. — and Lg — L,- cases. 

The Z' in gauged L^—Li {i = fx, r) cannot be sought for at the LHC since they only couple 
to leptons. However, they can be produced in e+e^ collisions and subsequently decay into 
e~^e^ or i^i^ pairs, and stringent contraints already exist from LEP/LEP2. The exchange 
of the Z' induces the following effective four-fermion interactions, relevant to e~^e~ colliders, 
among the charged leptons at energies far below the Z' mass: 
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The LEP collaborations fit their data to 

47r 
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(15) 



(16) 



with the 95% confidence limits on A± shown in Table [B The strongest constraint for the 
Le — case comes from the e+e^ fi^fi^ channel of L3, which translates to 

Mz' 



> 4.1 TeV 
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FIG. 2: Diagrams that contribute to neutrino oscillation matter effects in (a) the gauged X = 
B — aLf, — f3L^ — 'yLj- model, i = {e, ;U,t}, / = {n, d, e}, and (b) topcolor assisted technicolor, 
/ = {uL,UR,dL,dR,eL,eR}. 

while that for the — case comes from the e^e^ r^r^ channel of DELPHI, which 
translates to 



Mz ' 

9z' 



> 4.5 TeV . 



(18) 



Though these are the 95% confidence limits while that given in Eq. (HM is the 99% limit, 
it is clear that the bound on ^ will not lead to any improvement of already existing bounds 
from LEP/LEP2. 



B. Gauged B - {ale + (31^ + -fLr) 

In Refs. Is, 3, y, 3 , extensions of the SM gauge group to Sf/(3)c x SU{2)l x f/(l)y x f/( 1)^ 
with X = B — {aLe + (3L^ + '-)L-r) were considered. Again, the motivation was to explain the 
observed pattern of neutrino masses and mixings. The cases (a, 7) = (0,0,3), (3,0,0), 
and (0, |, |) were considered, respectively, in Refs. 3, [3], and In all cases, the condition 

a + /3 + 7 = 3 (19) 

is required for anomaly cancellation within the SM plus right-handed neutrinos^. When 
7^ 7^ 7) the U{l)x gauge boson, i.e. the Z' , couples to the three lepton generations 
differently, and can lead to extra neutrino oscillation matter effects. As in the gauged 
Le — Li case, the Higgs sectors of these models also necessarily distinguish among the lepton 
generations, but we relegate the discussion of their effects to section V. 



^ Only the right-handed neutrinos with non-zero X charge need to be included for anomaly cancellation. 
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For generic values of 7), the Z' couples to the quarks and leptons as 

Lz, = g^,r^Z'^, (20) 

where 

f 



3 

9 



(21) 



The forward scattering amplitude of the left-handed neutrino {i = e, fj,,T) on matter 
fermion F {F = p, n, e) due to Z'-exchange (c/. Fig. [2ti) is 

iMf = i+tgz'X^Ji+tgz') {ueLlm^i^iWiL) {j^^ {F\r^\F) . (22) 

Again, we can assume that the matter fermions are non-relativistic, so that only the time-like 
components of the currents need be considered. Then, we can make the replacements 

(e| Jx |e) = —a (e| e^e |e) — *• —aN^ , 

{V\ J'x \P) = I {P\ («^« + d^d) \p) i (2iVp + AT,) = AT, , 

{n\ Jx \n) = ^ (n| {u^u + dU) \n) ^ \{Nn + 2Nn) = Nr^ , (23) 

and 

{i^iL\l^a^i^iWeL) = M [^Il^il + ^m'^eR) WiL) = Mi'lL^iLWeL) , (24) 



which gives us 



iMf = -iX.,^^ {4 A.) (Xj^Np) , (25) 

Z 

where we have defined = X„ = 1. Summing over F = p,n, e, we find: 
iM = i 'Y 

F=p,n,e 

= -^^-.^(0U.J(^P + ^n-«Ar,) = -^K, (0lAj ' 

(26) 

where 



= +^-. it i^n + N,- aN,) (27) 
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can be identified as the effective potential experienced by the left-handed neutrino v^l as it 
travels through matter. Since the Earth is electrically neutral and is mostly composed of 
lighter elements, we can make the approximation A'^ ^ Np = = N , in which case 



(28) 



The effective ^ is then 



-4(a-2)(/5-7) 



{gz'/Mz'f 
{g/My^^ ■ 



(29) 



Vcc yyi^viw) 
When a = 2, the contribution of the matter electrons is cancelled by those of the matter 
nucleons and ^(2,/3,7) vanishes, regardless of the values of (3 and 7. When /3 = 7, the matter 
effects on i/^ and Vr will be the same, again resulting in C,(a,i3,/3) = 0, regardless of the value 
of a. 

In Fig. [3l we plot the dependence of C,z' on the Z' mass for selected values of gz' for the 
case a = /? = 0, 7 = 3, namely, the Z' couples to B — 3L^. In this case 



^(0,0,3) = -24 



{gz'/Mz 



6 



9z' 



(30) 



{g/Mw^ V2Gf \Mz' 

Ignoring the possible contribution of the Higgs sector, a bound on ^ of |^| < ^0 = 0.005 from 
Eq. (IHl) translates into: 

> ^ 8500GeV. (31) 



9z' 



V2Gf^o 



a=/3=0, 7=3 




2 3 4 
Mz' (TeV) 



FIG. 3: (,z' dependence on the Z' mass for the special case a = /3 = 0, 7 = 3. 
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FIG. 4: Lower bounds on Z' mass. 
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(0,0,3) 


0.65 


580 GeV 


~ 1 TeV 


5500 GeV 




0.35 


220 GeV 


~ 0.6 TeV 


3000 GeV 


V'-' 2' 2/ 


0.65 


500 GeV 
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0.35 
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TABLE II: Current and possible lower bounds on the Z' mass in gauged B — aL-^ — f3L^ — 
models. 

More generically, the bound on the Z' mass is 
Mz' 



> 



9z' 



yS^^^^ ^ v/|(«-2)(/3-7)|x(3500GeV) 



(32) 



This bound is plotted in Fig. [Has a function of (3 for three different values of gz', and two 
different values of a. The value of 7 is fixed by the anomaly cancellation condition, Eq. f|T9l) . 
to 7 = 3 — a; — /3. The region of the (/3, Mz') parameter space below each curve will be 
excluded. 



Let us now look at existing bounds. We limit our attention to the a = case, i.e. the 
Z' couples to i? — [3L^ — jLr, with /3 + 7 = 3. In this case, the Z' can be produced in pp 
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collisions and subsequently decay into yU+yU" or r+r~ pairs. The exchange of the Z' in this 
case leads to the following four-fermion interactions, relevant to colliders, between the 
charged leptons and the light quarks at energies way below the Z' mass: 

DO has searched for the contact interaction 

^ = ^T2 (^^'^^ + ^^''f^) (/i7/./w) (34) 



in its dimuon production data |27l] and has set a 95% confidence level limit of 

A+>6.88TeV. (35) 

This translates into 

^^'>V^x(l.lTeV). (36) 



9z' 



CDF has searched for the production of a Z' followed by its decay into t^t pairs 28] and 
has set a 95% confidence level lower bound of 

Mz> > 400 GeV (37) 

for a sequential Z' (i.e. a Z' with the exact same couplings to the fermions as the SM Z). 
Rescaling to account for the difference in couplings, we estimate 

^>V^x(lTeV). (38) 

9z' 

Limits on this model also exist from a global analysis of loop effects in LEP/SLD data j^, 
but they are weaker than the direct search limits from the Tevatron. In Table HTl we compare 
the bounds from LEP/SLD, CDF/DO, and the potential bounds from a measurement of ^ 
for two choices of (a, /?, 7), and two choices for the value of gz'- For the (a, P, 7) = (0, 0, 3) 
case, we can expect a significant improvement over current bounds. 

The sensitivity of the LHC to Z's has been analyzed assuming Z' decay into e~^e~ or 
/i"^yU~ pairs, or 2 jets 3^. For a sequential Z', the LHC is sensitive to masses as heavy 
as 5 TeV with 100 fb~^ of integrated luminosity. The Z' of the 7) = (0,0,3) model, 

however, decays mostly into r+r~, which will not provide as clean a signal as decays into 



the lighter charged lepton pairs. Ref. [lO| estimates that if gz' ~ (?' ~ 0.35, then the LHC 
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reach will be up to about 1 TeV with 100 fb~^. If this estimate is correct, the potential 
bound on Mz' from neutrino oscillation may be better than that from the LHC A complete 
detector analysis may show that the actual reach of the LHC is somewhat higher, but even 
then we can expect the neutrino oscillation bound to be competitive with the LHC bound 
for the (0, 0, 3) model. 



C. Topcolor Assisted Technicolor 

Another example of a model with a Z' which distinguishes among different generations 



is topcolor assisted technicolor 
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12[ |. Models of this class are hybrids of topcolor and 
technicolor: the topcolor interactions generate the large top-mass (and a fraction of the W 
and Z masses), while the technicolor interactions generate (the majority of ) the W and Z 
masses. The models include a Z' in the topcolor sector, the interactions of which helps the 
top to condense, but prevents the bottom from doing so also. To extract the interactions of 
this Z' relevant to our discussion, we need to look at the model in some detail. 

Though there are several different versions of topcolor assisted technicolor, we consider 
here the simplest in which the quarks and leptons transform under the gauge group 

^f/(3), X SU{?>)^ X U{l)s X U{1)^ X SU{2)l (39) 

with coupling constants g^si Qsw, du, Qiw, and g. It is assumed that g^g ^ g^^ and gu ^ giw 
SU{2)l is the usual weak-isospin gauge group of the SM with coupling constant g. The 
charge assignments of the three generation of ordinary fermions under these gauge groups 
are given in Table IIIII Note that each generation must transform non-trivially under only 
one of the S'f/(3)'s and one of the [/(l)'s, and that those charges are the same as that of 
the SM color, and hypercharge Y (normalized to Qem = -^3 + ^)- This ensures anomaly 
cancellation. 

At scale A ~ 1 TeV, technicolor, which is included in the model to generate the W and Z 
masses, is assumed to become strong and generate a condensate (of something which is left 
unspecified) which breaks the two S't/(3)'s and the two f/(l)'s to their diagonal subgroups: 

SU{?>)s X SU{?>)^ ^ SU{?>), , U{l)s X f/(l)^ U{1)y , (40) 

which we identify with the usual SM color and hypercharge groups. The massless unbroken 
SU(3) gauge bosons (the gluons G^) and the massive broken SU(3) gauge bosons (the so 
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TABLE III: Charge assignments of the ordinary fermions. The U{1) charges are equal to the SM 
hypercharges normahzed to Qem = l3 + Y- 

called colorons C^) are related to the original SU{3)s x SU{3)w gauge fields and X^^ 

by 



= Xsi^ COS 03 - Xy^ij, sin 
Gn = Xsn sin 63 + X^f, cos 03 



(41) 



where we have suppressed the color indices, and 



tan 03 



93u 
93s 



(42) 



The currents to which the gluons and colorons couple to are: 

93sJ3s^st, + 93wJ3w^wn = 93 (cot ^3 Jg^^ - tan 03J^J + g3 ( + J^J Gf, , (43) 

where 



11 1 

2 2~ ~^ 2 ■ 

93 93s 93w 

Since the quarks carry only one of the SU{3) charges, we can identify 
as the QCD color current, and g3 as the QCD coupling constant. 



(44) 



(45) 
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Similarly, the massless unbroken U(l) gauge boson and the massive broken U(l) gauge 
boson Z'^ are related to the original C/(l)s x C/(l)^ gauge fields Yg^i and by 

Z'^ = Ys^, cos Qx - Ywii sin Q\ 

= Y,^ sin ^1 + Y^^ cos (46) 



where 



tan^i = ^. (47) 



The currents to which the and Z'^ couple to are: 

guJ^Y,^ + giv^JLYn^f. = 91 (cot Jf, - tan^i Jf^) + g, ( Jf, + Jf^) 5^ , (48) 

where 

111 
91 9is 9tw 

Again, since the fermions carry only one of the U{1) charges, we can identify 

= jr, + j(; (50) 

as the SM hypcrcharge current, and gi as the SM hypercharge coupling constant g'. Note 
that the interactions of the colorons and the Z' with the third generation fermions are strong, 
while their interactions with the first and second generation fermions are weak. This results 
in the formation of a top-condensate which accounts for the large mass of the top quark. ^ 
Therefore, the interaction of the Z' in this model with the quarks and leptons is given by 

C^g' (cot^i Jf, - tan^i Jf^) Z'^ , (51) 

where g' is the SM hypercharge coupling, and 

J^s = I {hl^tL + h^bL) + ^tR-f^R - ^bR^bR - ^ {tlI^Tl + yrLl'^^rL) - Tr^Tr , 

1 2 11 

Jiyj = g {clI^cl + sli'^sl) + -crI^cr - -srj^sr - - {[IlYi^l + - ftRl^l^R 

1 - 2 1-1 

+ g {ulI^Ul + dL^dL) + -UrYur - -dR^dR - - {clI^Cl + ^eLl^l^eL) - eR^f'^CR . 

(52) 



^ The Z'-exchange interaction in the ti channel is attractive, but that in the bb channel is repulsive. This 
repulsion is assumed to be strong enough to counter the attraction due to the colorons and prevent the 
bottom from condensing. 
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FIG. 5: ^TT dependence on the Z' mass in the top color assisted technicolor model. 

The exchange of the Z' leads to the current-current interaction 

1 



(cot Q\J\s — tan^iJi^) (cot Q\J\s — tan^iJ; 



(53) 



the J\sJ\s part of which does not contribute to neutrino oscillations on the Earth, while the 
Jiwhw part is suppressed relative to the Jiwhs part by a factor of tan^^^i ^ 1. Therefore, 
we only need to consider the JuJiw interaction which only affects the propagation of 
{cf. Fig. [2]d). The forward scattering amplitude of against fermion F = p,n,e is given 
by 



-ig' cot 9i){+ig' tanOi) {u^lI ( i^rl^ Pli^t^ Wtl) 



x(F| 



,'2 



2M|, 



-Pl -PR]e 
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,/2 



2Ml, 
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n 



l(l + ^)(2iV, + iV„ 



4 4 4 



1 (\ 1 



2 V 6 



(iVp + 2iV„ 



1 



SMI, 



J2 X ^ 



Ml, 



) 



Note that the angle 6i has vanished from this expression and the only unknown parameter 
here is the Z' mass. 

The effective potentials felt by the different neutrino flavors are 



(54) 







N g 



.12 



8 Ml, ' 



(55) 
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and the effect ivie ^ is 

The dependence of ^tt on the Z' mass is shown in Fig. [51 The hmit \^tt\ < = 0.005 in 
this case translates to: 



sin^ 6i 



Mz' >Mz\l ^ ^ 440 GeV . (57) 

This potential limit from the measurement of ^ is much weaker than what is already available 
from precision electroweak data [l^, or from the direct search for pp —>■ Z'X t^t^X at 



CDF mentioned earlier 



28 



III. GENERATION NON-DIAGONAL LEPTOQUARKS 

Leptoquarks are particles carrying both baryon number B, and lepton number L. They 
occur in various extensions of the SM such as Grand Unification Theories (GUT's) or Ex- 
tended Technicolor (ETC). In GUT models, the quarks and leptons are placed in the same 
multiplet of the GUT group. The massive gauge bosons which correspond to the broken 
generators of the GUT group which change quarks into leptons, and vice versa, are vector 
leptoquarks. In ETC models, the technicolor interaction will bind the techniquarks and 
the technileptons into scalar or vector bound states. These leptoquark states couple to the 
ordinary quarks and leptons through ETC interactions. 

The interactions of leptoquarks with ordinary matter can be described in a model- 



independent fashion by an effective low-energy Lagrangian as discussed in Ref. 31] . Assum- 
ing the fermionic content of the SM, the most general dimensionless SU{3)c'xSU{2)l x U{1)y 
invariant couplings of scalar and vector leptoquarks satisfying baryon and lepton number 
conservation is given by: 

^ = ^F=2 + ^F=0 5 (58) 
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where 



F=2 



Si + giR 



Si 



+93L 



S, 



+ 



921 d'R-f^iL + 92Rqll^eii 



2fi + 92L 



V^2m + h.c. 



(59) 



h2LU]ih + h2RqEi'T2eR S2 + /12L 
+ hiLqlYh + hiRd^YeR 
+h3L \qlT7^h] % + h.c. . 



d 



R^L 



s. 



Vi^ + hiR 



URl^ej^ 



(60) 



Here, the scalar and vector leptoquark fields are denoted by S and V , respectively, their 
subscripts indicating the dimension of their SU{2)l representation. The same index is 
attached to their respective coupling constants, the g^s and /I's, with the extra subscript L 
or R indicating the chirality of the lepton involved in the interaction. For simplicity, color, 
weak isospin, and generation indices have been suppressed. The leptoquarks Si, Si, S3, V2, V2 
carry fermion number F = 3B + L = —2, while the leptoquarks 5*2, S2, Vi, Vi, V3 have F = 0. 

Rewriting the fermion doublets and the leptoquark multiplets in terms of the individual 
component fields, Eqs. (15^ and fp7]) are expanded as follows: 



-F=2 



9iLiuleL - dlu^) + gmiuR^R) S^ + gm 



^R^R 



+ 



g2L{d%-i^eL) + g2R{dli''ej^) 



^2^ + 



^2l(c?r7'''^l) + ^2i?«7%) 



2/i 



+g2L 
+g3L 



F=0 



h2L{uReL) + h2R{uLeR) 



st + 



h2L{uRVL) - h2R{dLeR) 



+ h.c. 



So 



(61) 



+h 



2L 



{dReL)S^ + {dRVL)S2 



+ 



hiiiuLl^VL + c^lt'^cl) + hiRidR-ff'eR] 



Vil + hlR 



UrY^r 



V2{uEreL)V^t + {uZl^VL - dL^^CL^l + V2{dLi^VL)V^} + h.c. .(62) 



Superscripts indicate the weak isospin of each field, not the electromagnetic charge. For 
fields with subscript 1, the superscript is redundant and may be dropped. The quantum 
numbers and couplings of the various leptoquarks fields are summarized in Table IIVI Note 
that the scalar 5*1 and the vector Vi^ do not couple to the neutrinos, so they are irrelevant 
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Leptoquark 


Spin 


F 


SU{3)c 


h 




Qem 


Allowed Couplings 


Si 







-2 


3 





1 
3 


1 
3 


gihitAer -dPrVr), giR(u%e r.) 




s9 





-2 


3 





4 
3 


4 
3 








1 


-2 


3 


2 
1 

2 


5 

6 


4 
3 

1 
3 


g2Lid'h^'^er), g2Rid'i^^er,) 
g2L(d^l^i^L)^ 92R(^l^eR) 




2/1 


1 


-2 


3 


+ k 

1 

2 


1 

6 


1 
3 

2 

Q 


g2L{V'%')^V^) 


S3 


53° 

S3 





-2 


3 


+1 

-1 


1 
3 


4 
3 

1 
3 

2 

3 


-V253L(^e^) 
-53L«ei + d5,z/i) 

\/2ff3L(^iiz^Z,) 




St 

'-'2 

So 
2 








3 


2 

_1 
2 


7 
6 


5 
3 

2 
3 


^2L (^^e L),h2R (uTeR ) 

h2L{uRl^L), —h2R{dLeR) 


S2 


St 

S2 








3 


2 

1 
2 


1 
6 


2 

3 

1 
3 


h2L{dReL) 

h2L{dRVL) 


Via 


V?, 


1 





3 





2 

3 


2 

3 


hiLiWl'^'^L + dL^^'eL), hiRidR'y^^eR) 






1 





3 





5 

Q 


5 

Q 


hlRiuRj^CR) 


V3^, 


V^3°. 


1 





3 


+1 

-1 


2 

3 


5 
3 

2 
3 

1 
3 


V^h3L(uL^^eL) 
h3L{uLl^l'L - dLl'^eL) 

V2h3L{dll^'UL) 



TABLE IV: Quantum numbers of scalar and vector leptoquarks with S'C/(3)c x SU{2)l x U{1)y 
invariant couplings to quark-lepton pairs (Qem = l3 + Y). 



to our discussion and will not be considered further. The isospin plus components of the 
remaining leptoquarks, namely S^, S^, S^, V^^, V^n, and V^, do not couple to the neutrinos 
either, but we will keep them in our Lagrangian since their coupling constants are common 
with the other components that do couple, and are important in understanding how the 
couplings are constrained by neutrinoless experiments. 

Since the leptoquarks must distinguish among different generation fermions to contribute 
to neutrino oscillation matter effects, we generalize their interactions by allowing the cou- 
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pling constants to depend on the generations of the quarks and leptons that couple to each 
leptoquark: 



F=2 



+ 



V, 



2fi 



+9iL 



u 



2/i 



F=0 



h2L{uiR^jL) - h^nidiLejR] 



+ h.c. 



(63) 



2L 



{diRejL)S^ + {diRUjL)S2 



+ 



+ h.c. . 
(64) 



Here, i is the quark generation number, and j is the lepton generation number. Summation 
over repeated indices is assumed. The interactions that contribute to neutrino oscillation 
matter effects are those with indices {ij) = (12) and (ij) = (13). It is often assumed in 
the literature that generation non-diagonal couplings are absent to account for the non- 
observation of flavor changing neutral currents and lepton flavor violation. However, the 
constraints from such rare processes are always on products of different (ij)- couplings and 
not on the individual non-diagonal couplings by themselves. For instance, non-observation 
of the decay efi constrains the product of (12) and (21) couplings, but not the 

(12) and (21) couplings separately, which allows one of them to be sizable if the other is 
small. Constraints on the individual (12) and (13) couplings actually come from precision 
measurements of flavor conserving processes, such as -Rtt = r(7r ^v^)/T{tx ez/g) which 
constrains the square of the (12) coupling, and those constraints are not yet that strong 32|. 

In the following, we calculate the effective value of ^ induced by the exchange of these 
leptoquarks. The leptoquark fields are naturally grouped into pairs from the way they couple 
to the quarks and leptons: {Si^S^), {82,82), (V2,V2), and (V^i,V^). We treat each of these 
pairs in turn, and then discuss the potential bounds on the leptoquark couplings and masses. 
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A. Si and S3 lepto quarks 

cO 




p + k 




Uj{k) 




S, 



(a) Q = 1, 3 



p + k 




v,{k) 



FIG. 6: Diagrams contributing to neutrino oscillation matter effects from the exchange of (a) 
Si or the isospin component of 5*3 , and (b) the isospin —1 component of S3. The EM charge 
Qem = H + Y for S^ and are +g, while that for S^ is — |. 

The {ij) = (12) and (13) interactions of the leptoquarks 5*1 and S3 are, respectively, 

^ = -9rLKi^,L)Si - 9llKiy,L)Si + h.c. , (65) 



and 



^ = 93L Hdlu,L)Sl^ + V2{ulu,L)S., 



+ 9ll 



h.c. 



(66) 



The interactions described by Eqs. fl65l) and fl66|) can be written in a common general form 

as 

£ = A {Y'Pl + \*{VPr q')S , (67) 

where q = u or d. The Feynman diagrams contributing to neutrino oscillation matter 
effects are shown in Fig. [61 At momenta much smaller than the mass of the leptoquark, the 
corresponding matrix element is 



(68) 



Using the Fiertz rearrangement 



{l7Pnq'){q-PLu) = -^{vYPLiy){q'l,PRq'') = +\ {V^ Pl i^) {qi ,Pl q) , (69) 



we obtain 



s s 
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where 



V= ^^'^ 



4 M| 



(71) 



Applying this expression to the Si case, the effective potential for the neutrino of generation 
number j is: 



The effective ^ is then 



Nd\g\i\ 

4 



{N^ + 2N^)\g\i 
4 



m2 



4 ' 



Csi — 



For the case, the effective potential is 



4 M2 2 ML 



SAT 



2Ml 



SI 



and the effective ^ is 



+3 



cc 



1 2 

+ 



w 



Ml, Ml_ 



(72) 



(73) 



(74) 



(75) 



In the case of degenerate mass, Mgo — Mg- = Mg^, we have 



>53 



+9 



(76) 



B. S2 and S2 leptoquarks 



The relevant interactions are 



for S2 and 



jC = h2L{uRU^L)S2 + h2L{^RVrL)S2 + h.C. (77) 



C = hll{dRU^L)S2 + /i2i(c?i?'^rL)52- + h.C. (78) 
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+i hli 



2L 



U{p) Vj{k) 



2L d{p) 
t * 



^2" i 



k — p 



k — p 



u{p) 



+i hli 



2L 



uj{k) dip) 



+i h\i 



2L 



(a) 



{b) 



FIG. 7: Diagrams contributing to neutrino oscillation matter effects from the exchange of (a) , 
and (b) §2 ■ The EM charge Qem = ^3 + 1" for 5*2" +§> '^hile that for §2 is — |. 



for 5*2 leptoquarks. Both fl77|) and flTHj) can be written in a common general form as 

C = \{qPLu)S + \*iVPRq)S , 



(79) 



where q = u 01 d. The Feynman diagram contributing to neutrino oscillation matter effects 
is shown in Fig. [7^. For momenta much smaller than the mass of the leptoquark, the 
corresponding matrix element is 



^M = {-lYW {y.q\ {VPRq) ^-^J (qPLiy) W,q) 
Using the Fiertz identity given in Eq. (I69l) again, we obtain 



(80) 



= -t-^ {u\u^^Plu\u) {q\q^,PRq\q) ^ -^^N^ {<Pl<Pu) = "^K (0t0j , (81) 
s s 



where 



4 M| 



(82) 



Applying this expression to the 5*2 case, the effective potential for the neutrino of generation 
number j is 



and the effective ^ is 



Nu\h'l\ 
' 4 ML 



{2Np + Nn) |/i2i| 



Mi 



3N\h'l\ 



4 M' 



2 ' 



(83) 



K3 - 



Vcc 



-3 



\hllf-\hlif)/Ml. 



(84) 



w 
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The effective potential for the 5*2 case is 



4 ML 4 ML 



4 M 



2 ' 



and the effective ^ is 



K3 - 
Vcc 



-3 



\m\')/Mi 



w 



(85) 



(86) 



C. V2 and ^2 

Uj{h) 




Vj{k) 




(a) 



(ft) 



FIG. 8: Diagrams contributing to neutrino oscillation matter effects from the exchange of (a) , 



and (b) ■ The EM charge Qem = h + Y for KT is +|, while that for K 



is 



'3- 



The relevant interactions for Vo are 



and those for are 



Both (1871) and (IHHl) can be written in a common general form as 



(87) 



(88) 



(89) 



The Feynman diagrams contributing to neutrino oscillation matter effects are shown in 
Fig. [HI For momenta much smaller than the mass of the leptoquark the corresponding 
matrix element is 



tM = {-iY\\\'{p,q\iV^^PLq') 



1 



(90) 
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Using the Fiertz rearrangement 



{vYPLq''){q'l^.PLiy) = {vrPLy){q''lt.PLq') = - {V^PlI^) {ll^^PRq) , (91) 
we obtain 

(H^7'^^L^k) {q\q^,PRq\q)^^^N,{<pl<p,)=-tKi<pi<p,) , (92) 



where 



N„\X 



2 



^ ■ (93) 



2 



Applying this to the case, the effective potential for the neutrino of generation number 
j is 



" 2 2 M2 ~ 2 • ^ ^ 

''2 ''2 

The effective ^ is 



The effective potential for the V2 case is 



The effective ^ is 



" 2 M.?_ ~ 2 M.?_ ~ 2 M? ■ ^ ' 



Vcc ' " gVM^ 



ever - ^^^V^ - +6 • " • (97) 



D. Vi and leptoquarks 



The relevant interactions for Vi are 



(98) 



and those for V3 are 



hll 



{uZ^'VrLWl + V2(dZ7'^i.,i)y3; 



+ /i.e. 



(99) 
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u{p) 



dip) 



p — k 



p — k 



u{p) 



aL 



Uj{k) 



dip) 



(a) Q = 1, 3 



ib) 



FIG. 9: Diagrams contributing to neutrino oscillation matter effects from the exchange of (a) 
or the isospin component of V^, and (b) the isospin —1 component of V^- The EM charges 
Qem = H + Y for Vi and are +|, while that for is — |. 

The interactions described by Eqs. fl98l) and fl99|) can be written in a common general form 



as 



C = \ [q^^^PL v)V + \*iVY'PL q)V . 



(100) 



The Feynman diagrams contributing to neutrino oscillation matter effects are shown in 
Fig. O For momenta much smaller than the mass of the leptoquark the corresponding 
matrix element is 



tM = i-t)'\X\'{iy,q\ivrPLq) 



M2 ' ^'^ 



iQlt^Pi^) W, q) 



Using the Fiertz identity given in Eq. (I9T1) again, we find 



where 



" ~ 2 MJ ■ 



Applying this result to the Vi case, effective potential is 



The effective ^ is 



+ 



I'^IL 



2 (M, 



Vi, 



(2iVp + iV„.) \h\i\ 
2 (Ml 



Vi, 



3N 

'^iMv,)' 



IVi 



(101) 



(103) 



(104) 



(105) 
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The effective potential for the V3 case is 



The effective ^ is 

^^3 



2 M^o 



M2. 



(2iVp + iV„) , (7Vp + 2iV„) 



2M2^o 



+ 



''3L 



+ 



^3° 



3Ll 



1 



In the case of degenerate mass, Myo = My- = My^, we have 



M,2 M?, 



V,- 



(106) 



(107) 



V3 



-18 



(108) 



All couplings = 0.5 




1 2 3 

Mlq (TeV) 



FIG. 10: ^LQ dependence on the leptoquark mass for ^AA|q = 0.5. (a) ^i; (b) V2, V2', (c) S3; 
(d) S2, S2; (e) Fi; (f) ^3- 
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LQ 


Clq 




upper bound from \^\ < 


current bounds from Ref. [321 


Si 


+3 


blip - blip 


1.1 X 10-3 


(5|f)2< 0.008 (R^) 
(5|i)2<0.7 (r^TTi/) 


S3 


+9 


bsil'-biP 

1 -'O-L/ 1 1 O-L/ 1 


3.7 X 10""^ 


{g\l? < 0.008 {R^) 
(<?^i)'<0.7 (r^vri.) 


S2 


-3 


11,12 |2 _ |;^13 12 
I"'2LI \^2L\ 


1.1 X 10^3 




S2 


-3 


|I12|2 _ |I13|2 


1.1 X 10^3 


{hllf<2 ifiN^f^X) 


V2 


+6 


b2il'-b2'iP 


5.5 X 10-"^ 


(<7^i)2<l ifiN^fiX) 


V2 


+6 


\~&-\& 


5.5 X 10""^ 


Cgliy<5 i^,N^^,X) 


Vi 


-6 


l^lil'-bliP 

1 ±i/ 1 1 11^ 1 


5.5 X 10-"^ 


{h\lf < 0.004 (R^) 
(/i|i)2<0.1 (D^f^u) 


% 


-18 


1^,12 |2 _ |;^13 12 


1.8 X 10-"^ 


{hllf<omA (R^) 

(/lli)2<0.1 (D^M^^) 



TABLE V: Constraints on the leptoquark couplings with all the leptoquark masses set to 100 GeV. 
To obtain the bounds for a different leptoquark mass Mlq, simply rescale these numbers with the 
factor {Mlq /WO GeV) 2. 

E. Constraints on the Leptoquark Couplings and Masses 



Assuming a common mass for leptoquarks in the same SU{2)l weak-isospin multiplet, 
the effective ^ due to the exchange of any particular type of leptoquark can be written in 
the form 

SXlo/Mln Clo /^^^io^ . ^ 



Here, C^g is a constant prefactor, and 5\\q represents 



2 _ I \ 12 |2 I \13 |2 



(110) 



where A^g is a generic coupling constant. The values of Clq and 5X\q for the different types 
of leptoquark are listed in the first two columns of Table |Vl In Fig. [10], we show how E,lq 



depend on the leptoquark mass Mlq for the choice ^ 5\\q = 0.5, where we have assumed 
SXlq > 0. To obtain the picture for the case when ^Aig < 0, the vertical axis of the graph 



27 



^ 4000 

o 

^3000 




FIG. 11: Lower bounds on the leptorquark masses, (a) Si, 5*2, 5*2; (b) Vi, V2, V2; (c) 5*3; (d) V3 



Process 




LQ 


Assumptions 


95% CL bound 


Reference 


pp ^ LQLQ X ^ {ju){ju)X 


(**) 


5 


p = o(") 


117 GeV 


CDF [34] 


pp^LQLQX ^{ju){jv)X 


(**) 


5* 


p = o 


135 GeV 


DO [35j 


pp^LQLQX ^{jf,){jf,)X 
pp^LQLQX ^{jfi){ju)X 


(*2) 


5* 


(3 = 0.5 


208 GeV 


CDF [36] 


pp^LQLQX ^{jfi){jfi)X 
pp^LQLQX ^{jfi){jiy)X 


(*2) 


S 


(3 = 0.5 


204 GeV 


DO [37] 


pp^ LQnX ^ {in)^ X 


(*2) 


S 


13 = 0.5, A = l('') 


226 GeV^'^) 


DO [38] 


pp ^ LQLQ X ^ {jT){jT)X 


(*3) 


V 


minimal coupling [40] 


251 GeV 


CDF [39] 



TABLE VI: Direct search limits on the Leptoquark mass from the Tevatron. is the assumed 
branching fraction B{LQ q£) = 1 — B{LQ qu), and ^^^^A is the Yukawa coupling of the 
Leptoquark with the quark- lepton pair. Combined bound with the pair production data. 



should be flipped. The constraint \^lq\ < ^0 translates into: 



Mlq > M, 



w 



\C 



LQ\ 



6 



\Clq\\6XIq\ 



CiQll-^AiQl X (1700GeV) . (Ill) 



The resulting bounds are shown in Fig. [TTl where the regions of the {Mlq, y |5A|q]) param- 
eter space below each of the lines will be excluded. One can also fix the leptoquark mass 
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and obtain upper bounds on the leptoquark couplings: 



The values when M^q = 100 GeV are listed in the third column of Table El The bounds 
for a different choice of leptoquark mass Mlq can be obtained by multiplying by a factor 
of {Mlq /lOOGeVy . This result can be compared with various indirect bounds from rare 
processes which are listed in the last column of Table |Vl As can be seen, the limits from 
1^1 ^ ^0 can significantly improve existing bounds. 

Limits on leptoquark masses from direct searches at the Tevatron are listed in Table IVIi 
Bounds from LEP and LEP II are weaker due to their smaller center of mass energies. 
Since neutrino oscillation is only sensitive to leptoquarks with (ij) = (12) and/or (ij) = 
(13) couplings, we only quote limits which apply to leptoquarks with only those particular 
couplings, that is, leptoquarks that decay into a first generation quark, and either a second 
or third generation lepton. Though it is usually stated in collider analyses that leptoquarks 
are assumed to decay into a quark-lepton pair of one particular generation, it is often the 
case that the jets coming from the quarks are not flavor tagged. Analyses that look for 
the leptoquark in the quark-neutrino decay channel are of course blind to the flavor of the 
neutrino. Therefore, the bounds listed apply to leptoquarks with generation non-diagonal 
couplings also. 

As can be seen from Table IVH the mass bounds from the Tevatron are typically around 
200 GeV and are mostly independent of the leptoquark-quark-lepton coupling A. This in- 
dependence is due to the dominance of the strong interaction processes, qq annihilation 
and gluon fusion, in the leptoquark pair-production cross sections, and the fact that heavy 
leptoquarks decay without a displaced vertex even for very small values of A: the decay 
widths of scalar and vector leptoquarks with leptoquark-quark-lepton coupling A are given 
by X'^Mlq/IQtt and X'^Mlq/^Ah, respectively, which correspond to lifetimes of O(10~^^) sec- 
onds for Mlq = O(IO^) GeV, and A = O(10~^). In contrast, the potential bound on Mlq 



from neutrino oscillation, Eq. flllll) . depends on the coupling a/IC'^qI |5A|q|, but can be 



expected to be stronger than the existing ones for y |Clq||5A|q| as small as 0.1. 

Bounds on leptoquarks with (ij) = (12) couphngs can also be obtained from bounds on 



29 



contact interactions of the form 

An 

C = (qYPxq) (/i7M^L/i) , (113) 

where X = L or R, and q = u or d. For instance, at energies much lower than the leptoquark 
mass, the exchange of the Si leptoquark leads to the interaction 

|ol2|2 

= ^2Mf (^^'^^^) if^l.PLfi) . (114) 

The remaining cases are listed in Table IVIII The 95% CL lower bounds on the A^'s from 
CDF can be found in Ref. ^^], and the cases relevant to our discussion are listed in Ta- 
ble IVIllI These bounds translate into bounds on the leptoquark masses and couplings listed 
in Table IVnl Clearly, the potential bounds from |^| < ^q, also listed in Table IVlTl are much 
stronger. It should be noted, though, that the results of Ref. |[26|] are from Tevatron Run I, 
and we can expect the Run II results to improve these bounds. Indeed, Ref. 27|] from DO, 
which we cited earlier in the Z' section, analyzes the Run II data for contact interactions of 
the form 

An - 

C = {uY'Pxu + d^^Pxd) (/i7M^L/i) , X = L or i? , (115) 

and places 95% CL lower bounds on the A^'s in the 4 ~ 7 TeV range. While these are 
not exactly the interactions induced by leptoquarks, we can nevertheless expect that the 
bounds on the A^^'s will be in a similar range, and thereby conclude that the Run II data 
will roughly double the lower bounds from Run I. Even then. Table IVIII indicates that the 
potential bounds from |^| < will be much stronger. 



The prospects for leptoquark discovery at the LHC are discussed in Refs. [30|, |4l|. At 
the LHC, leptoquarks can be pair-produced via gluon fusion and quark-antiquark anni- 
hilation, or singly-produced with an accompanying lepton via quark-gluon fusion. The 
pair-production cross section is dominated by gluon fusion, which does not involve the 
leptoquark-quark-lepton coupling A, and is therefore independent of the details assumed for 
the leptoquark interactions. Once produced, each leptoquark will decay into a lepton plus 
jet, regardless of whether the coupling is generation diagonal or not. The leptoquark width 
in this decay depends on A, but it is too narrow compared to the calorimeter resolution 
for the A-dependence to be of relevance in the analyses. Therefore, though the analyses of 



Refs. 



30 



4l| assume specific values of A and generation diagonal couplings, we expect their 



conclusions to apply equally well to different A-values and generation non-diagonal cases: 
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LQ 



5i 



S2 



5, 



Vi 



V2 



Induced Interaction 



5*2 



+ 



MP 



TUT 



+ 



V2 



T/jTsp : 



CDF 95% CL [26] 



Ms, 
\9ll\ 



> 0.68 TeV 



Ms, 



> 0.72 TeV 



M, 



52 



> 0.38 TeV 



My, 

I/) 12 I 



> 0.48 TeV 



My, 
\9ll\ 



> 0.56 TeV 



V2 



> 0.85 TeV 



> 3.0 TeV 



IL 



Ms, 



> 3.0 TeV 



M. 



^ > 3.0 TeV 



> 5.2 TeV 



My, 



> 4.3 TeV 



My, 



> 4.3 TeV 



2L 



> 4.3 TeV 



2L 



> 7.4 TeV 



TABLE VII: The quark-muon interactions induced by leptoquark exchange, and the bounds from 
CDF (26I compared with potential bounds from neutrino oscillations. Only the couplings that also 
contribute to neutrino oscillation are listed. Analysis of the Tevatron Run II data is expected to 
improve the CDF bound by a factor of two. 



(qn) chirality 


^tf. (TeV) 


A^^ (TeV) 


(TeV) 


A,^ (TeV) 


(LL) 


3.4 


4.1 


2.3 


1.7 


(RL) 


3.0 


3.6 


2.0 


1.9 



TABL 
DO 



VIII: The 95% CL lower bound on the compositeness scale from CDF 



Results from 



27| | do not provide limits for cases where the muons couple to only u or d, but we expect the 



bounds to be in the range 4 ~ 7 TeV. 



for (3 = B{LQ qi) = 0.5, the expected sensitivity is up to M^q ^ 1 TeV with 30 ^ fb of 
data [41]. Again, in contrast, the the potential bound from neutrino oscillation, Eq. (lllip . 



depends on the coupling y IClq] 1(5A|q1. If y IClqUSX^^qI = 0(1), then Eq. (11 111) will be 
competitive with the expected LHC bound. 
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IV. SUSY STANDARD MODEL WITH R-PARITY VIOLATION 



Let us next consider contributions from R-parity violating couplings. Assuming the 
particle content of the Minimal Supersymmetric Standard Model (MSSM), the most general 
R-parity violating superpotential (involving only tri-linear couplings) has the form 13 1 



Wft = ^\ijkL,L,Ek + KjkUQjDk + \K,kUiDjDk , (116) 



where Lj, Ei, Qi, Di, and Ui are the left-handed MSSM superfields defined in the usual 
fashion, and the subscripts /c = 1,2,3 are the generation indices. (Note, however, that 
in some references, such as Ref. [14|, the isospin singlet superfields Ei, Di, and Ui are 
defined to be right-handed, so the corresponding left-handed fields in Eq. ( 11161) appear 
with a superscript c indicating charge-conjugation.) SU(2)l gauge invariance requires the 
couplings \ijk to be antisymmetric in the first two indices: 

whereas SU (3) gauge invariance requires the couplings X'-jj. to be antisymmetric in the latter 
two: 

Kjk = ~Kkj ■ (118) 
These conditions reduce the number of R-parity violating couplings in Eq. (11161) to 45 (9 
Ajjfc, 27 A-j^, and 9 A^'^-;,). The purely baryonic operator UiDjDk is irrelevant to our discussion 
on neutrino oscillation so we will not consider the A'/^-^ couplings further. We also neglect 
possible bilinear R-parity violating couplings which have the effect of mixing the neutrinos 
with the neutral higgsino: their effect on neutrino oscillation has been discussed extensively 



by many authors IJ, |42|, |43 1 . 



A. LLE couplings 

The LLE part of the R-parity violating Lagrangian, Eq. (I116p . expressed in terms of the 
component fields is 

Clle = Xijk[i>iLe^ejL + Cjie^^^iyiL + eljiU-L(^jL] + h.c. (119) 

The second and third terms of this Lagrangian, together with their hermitian conjugates, 
contribute to neutrino oscillation matter effects. The corresponding Feynman diagrams are 
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FIG. 12: LLE interactions that contribute to neutrino oscillation matter effects.. 



shown in Fig [121 Since Xijk is antisymmetric under i ^ j, it follows that i ^ j- Calculations 
similar to those for the scalar leptoquarks yield 



iVe 
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M? ^ 



M? 



(120) 



or if we write everything out explicitly: 
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(121) 



The effective ^ is 
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|A2jl 



E |A3jl| \ ^ 



|A2ij|^ - lAaijl^ 



w 



(|A 



j=l,3 ejL j=l,2 ^jL 
eiii 
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iA3iir) 



+ |A 
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Ml 
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IA212I ^ I A 
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Ml 
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M? 
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(122) 



For degenerate s-electron masses Mg^.^ = Mg^.^ = Mg^. , we have 



6 



|A 



2311 



IA122I ~ |A 



1321 



IA231I IA123I + I A 



133 I 



where we have used Xijk 



M? 

62 



Ml 



—Xjik to reorder the indices. 



(123) 
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FIG. 13: LQD interactions that contribute to neutrino oscillation matter effects.. 
B. LQD couplings 

The LQD part of the R-parity violating Lagrangian expressed in terms of the component 
fields is 



'LQD 



X 



ijk 



(iihdkRUjL + UjidkuCiL + dlj^e'^i^UjL 



+ h.c. 



(124) 



The second and third terms of this Lagrangian, together with their hermitian conjugates, 
contribute to neutrino oscillation matter effects. The corresponding Feynman diagrams are 
shown in Fig [131 Calculations similar to those for the scalar leptoquarks lead to the following 
effective potential for neutrino flavor uf. 



l^^lil 



The effective ^ is 



M. 



^3N f\X',,,\' 



A' 



4 \ M, 



r/2 



M, 



Vcc 



-3E 



I ^31i I 



/M? 



For degenerate d-squark masses M^,^ = Mg,^ = M^, , we have 



(125) 



(126) 



A il^2il| - |^3il| + |^21i 



I A' 



31j\ 



/M; 



(127) 
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FIG. 14: Dependence of on the smuon, sdown, and masses for y*^-^? = 0-5 in the (a) 
LLE R-parity violating interaction; (b) LQD R-parity violating interaction; and (c) the Zee/Babu- 
Zee models. 

C. Constraints on the R-parity Violating Couplings and Squark/Slepton Masses 

To illustrate our result for R-parity violating interactions, we simplify the analysis by 
assuming that only the A122 and A132 couplings are non-zero for the LLE case, and only the 
A211 and X'^ii couplings are non-zero for the LQD case. Under these assumptions, only the 
smuon, 62 = ft, contributes in the first case, and only the sdown, di = d, contributes in the 
latter. The corresponding ^'s are 



Cm 



+ 



where 





- 1 ( 




J a' a 


Ay/2GF \ 

6 


-.Mi) 


{g/MwY ~ 




V'l 


^^fi = -^122 P 


~ 1^1321 5 




= -^211 P 







128) 



(129) 

Fig. [13] shows how and depend on masses of the smuon and the sdown for a specific 
choice of couplings: ^(5A? = ^(5A| = 0.5 (we have assumed 5A| and 5A? to be positive). 
The bound |.^| < .^o = 0.005 translates into: 




5\l\ X (1700 GeV) 
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FIG. 15: Lower bounds on (a) the smuon mass in the LLE R-parity violating interaction model, 
(b) the sdown mass in the LQD R-parity violating interaction model, and (c) the mass in the 
Zee/Babu-Zee models, respectively. 




6 



X (4300 GeV) 



(130) 



The resulting graphs for the lower mass bounds are shown in Fig. [151 The regions of the 



M, 



ill 



and (Mj, 



(5Aj| ) parameter spaces below each of the lines are excluded. 



One can also fix the smuon and sdown masses and obtain upper bounds on the R-parity 
violating couplings: 



< 



'4:V2GFioM-^ = (0.057) 



Mr, 



<5A|| < 



'4^2^^ 
6 



(0.023) 



100 Gev) 



100 GeV 



;i31) 



These relations are actually more useful than Eq. (11301) since if the smuon and sdown exist, 
their masses will be measured/constrained by searches for their pair-production at the LHC, 
independently of the size of possible R-parity violating couplings. 

Current bounds on R-parity violating couplings come from a variety of sources HQ. 
The current indirect bounds of the four couplings under consideration from low-energy ex- 
periments are listed in Table IIX[ Comparison with Eq. (11311) shows that the bounds on 
Ai22 and A 132 are already fairly tight, and neutrino oscillation will do little to improve them. 
On the other hand, the bounds on A211 and can potentially be improved by factors of 
roughly 2.5 and 5, respectively. 
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Coupling 



n32| 



l'^122'^132l 



2111 



I A', 



3111 



Current 2a Bound 



0.05 



100 GeV 



0.07 



M, 



100 GeV 



(2.2 X 10" 



100 GeV 



0.06 



100 GeV 



0.12 



100 GeV 



Observable /Process 



Vud from nuclear (3 decay /muon decay 



r(r e i/gZ^r) 



r(T- 



3^ 



r vr- 



e Ue. 



T{tt- 



T(7r- 



TABLE IX: Current 2a bounds on R-parity violating couplings from Ref. ^14] . These bounds 
assume that each coupling is non-zero only one at a time. 

Bounds on R-parity violating couplings from ep and pp colliders come from searches 
for s-channel resonant production of sparticles. The bounds from the ep collider HERA 
necessarily involve the couplings Aj ■ since^the squark must couple to the first generation 



lepton (electron or positron) |44l . I45l . I46l . I47j so we will not discuss them here. The bound 
from the Tevatron comes from the analysis of DO which looked for the R-parity violating 
processes du ^ ft or dd ^ z>^, which occur if ^ 0, followed by the decay of the slepton 
via the R-parity conserving processes /2 ^ Xi 2 3 4 '^a* ^ Xi^2 A'- The neutralinos and 
charginos produced in these processes cascade decay down to the Xi (the assumed lightest 
supersymmetric particle, or LSP) which decays via a virtual smuon, muon-sneutrino, or 
squark though the R-parity violating \2ii coupling again into a muon and two jets, giving 
2 muons in the final state. The bound on the value of from this analysis depends in a 
complicated manner on all the masses of the particles involved in the processes. If one uses 
a minimal supergravity (mSUGRA) framework 49| with tan/3 = 5, /i < 0, and Aq = 0, then 
the 95% bound is X211 < 0.1 assuming = 363 GeV [48^. A similar bound would result 
from Eq. f|X31l) if = 460 GeV. However, since squarks are generically much heavier than 
sleptons l^toj], the existing DO bound is effectively stronger than the potential bound from 
lel < Co. 
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FIG. 16: Diagrams which generate the Majorana masses and mixings of the neutrino in the (a) 



Zee 



50l ] and (b) Babu-Zee 



5l[ | models. 



V. EXTENDED HIGGS MODELS 



Most models, including the Standard Model (SM) and its various extensions, possess 
Higgs sectors which distinguish among the different generation fermions. The models dis- 



cussed in section |TT] are necessarily so, and so are the Zee 50(] and Babu-Zee 5]| models 
of neutrino mass, as well as various triplet Higgs models [52|]. As representative cases, we 
consider the effect of the singlet Higgs in the Zee and Babu-Zee models, and that of a triplet 
Higgs with hypercharge Y = +1 {Qem = -^3 + 



A. Singlet Higgs in the Zee and Babu-Zee Models 



In the Zee 50j and Babu-Zee 51j models, an isosinglet scalar with hypercharge y = +1 
is introduced, which couples to left-handed lepton doublets as 



A = Xab{ilLCi(T2ibL)h+ + h.c. = Xab {ilLia2 ibl) + h-C. 



(132) 



where {ab) are ffavor indices: a,b = e, /i, r. The hypercharge assignment prohibits the 
fields from having a similar interaction with the quarks. Due to SU{2) gauge invariance, the 
couplings Xab are antisymmetric: Xab = —Xba- This interaction is analogous to the R-parity 



violating LLE coup 
In the Zee model 



ing with playing the role of the slepton. 



50|, in addition to the /i^, two or more SU{2) doublets 0^ (a = 1, 2, ■ ■ ■ ) 



with hypercharge Y = ~\ are introduced which couple to the via 

C^cj,h = Map {(pl iT2 (f)p)h^ + h.c. , 



(133) 
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and to the fermions in the usual fashion. The couphngs Map are antisymmetric, just hke 
Xab, which necessitates the introduction of more than one doublet. In this model, Majorana 
masses and mixings of the neutrinos are generated at one-loop as shown in Fig. [TBk . The 
extra doublets can also contribute to neutrino oscillation depending on their Yukawa cou- 
plings to the leptons, but we will assume that their effect is negligible compared to that of 
the h^. 

In the Babu-Zee model 511, in addition to the h^, another isosinglet scalar k~^^ with 
hypercharge Y = +2 is introduced which couples to the right-handed leptons and via 

= X'ab -Mh^h+k- + h.c. , (134) 

where A^^ = A'^^^^. In this model, Majorana masses and mixings of the neutrinos are generated 
at the two-loop level as shown in Fig. [T6b . In this case, the extra scalar, k, does not 
contribute to neutrino oscillation. 
Expanding Eq. (11321) . we obtain 

^ = 2 [ Ae^ ( t^/^L - i^e^ ) + Aer ( I^T^ " KlC^ ) + A^^ ( I^T^ " I^/^L ) ] + ^-C. 

(135) 

Keeping only the terms that are relevant for neutrino oscillation matter effects, we have 

- 2 ( Ae;.I^e^ + XerK^ej^)h+ + h.c. (136) 

The corresponding Feynman diagram is shown in Fig. [T71 
Calculations similar to those for the Si leptoquark yield 

and 

, _ - _ AM-\Xer\')/Ml _ 1 f6Xl\ 

" Vac {g/MwY ^V2G^ Wh) ' ^ ^ 

where we have defined 6X1 = l-^e/^P ~ |AerP- The dependence of C,h on the h"^ mass is 
plotted in Fig. [Tl]for the case ■>/ SX\ = 0.5, where we have assumed 6X1 > 0- The bound 
1^1 < ^0 = 0.005 translates into 

6X1 



Ml 



< V2Gf^o = (8.2 X 10"^) GeV-^ (139) 
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FIG. 17: Contribution to neutrino oscillation matter effects from a singly-charged Higgs in the 
Zee, Babu-Zee, and Y = 1 Triplet Higgs models. 



or 



\S\l\ X (3500 GeV) . 



(140) 



V2Gf ^0 

This result is represented graphically in Fig. [151 The region of the (M/i, a/|5A^|) parameter 
space below the constructed line would be excluded. 

A constraint on the exact same combination of the couplings and mass of the as 
above exists from r decay data: The measured value of the r~ — Ure'Ue branching fraction 
imposes the constraint [531] 



Ml 



< (3.4 X 10"^) GeV"' , 



;i41) 



which is clearly stronger than Eq. (11391) . 



B. Triplet Higgs with Y = +1 

We denote the components of an isotriplet Higgs with hypercharge F = +1 as 

A+ 
AO 



(142) 



It is customary to write this in 2 x 2 matrix form: 



A" 



V2 



+ A+fTs + A^ 



cTi + ia2 



A+/V2 A++ 
AO -A+/V2 



(143) 



The coupling of A to the leptons is then 

£a = V2K,{f,^Cta2Ai,^)+h.c. = {I^^ta^ Ai,^) + h.c. (144) 
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Model 


Stronger than existing bounds? 


Competitive with LHC? 


Gauged — and Lg — Lj. 


No 




Gauged B — 3Lt- 


Yes 


Yes 


Top color Assisted Technicolor 


No 




Lepto quarks 


Yes 


Yes* 


R-parity violation 


No 




Zee, Babu-Zee, Triplet Higgs 


No 





TABLE X: The result of our survey. The potential bound from |^| < .^o = 0.005 is compared with 
existing bounds, and the expected bounds from the LHC. If the existing bound is already stronger, 
no comparison with the LHC bound is made. *The leptoquark bound will be competitive with the 



LHC, provided that ^J\CLQ\\SXlQ\ = 0(1). 

This time, the couphngs are symmetric in the flavor indices A'^^ = A'^,^, and the factor of 
is thrown in for latter convenience. Expanding out, we find 



A 



ab 



V2 {yl^v,^) A° - [vl^e,^ + e^^j/,^) A+ - (e^^e,^) A++ 1 + h.c. (145) 



and the terms relevant to neutrino oscillation in matter are: 



- 2 ( AL vl^ + A',^ + A',, u';^ ) A+ + h.c. 



(146) 



Of these, the Agg term does not affect ^, while the other terms are precisely the same as those 
listed in Eq. fll36p . So without further calculations, we can conclude that all the results of 
the previous subsection apply in this case also. 



VI. SUMMARY AND CONCLUSIONS 

In this paper, we surveyed the potential constraints on various models of new physics 
which could be obtained from a hypothetical Fermilab^HyperKamiokande, or similar type 
of experiment. We assumed that the parameter ^, defined in Eq. (jlj), could be constrained 
to 1^1 < ^0 = 0.005 at the 99% confidence level. This places a constraint on the couplings 
and masses of new particles that are exchanged between the neutrinos and matter fermions. 

Table 1X1 summarizes our result. Of the models surveyed, the potential bound on gauged 
B — SLj. from |^| < can be expected to be stronger than the expected bound from the LHC. 
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Bounds on generation non-diagonal leptoquarks can be competitive if y |C'lq||(^A|q| = 0(1). 
For these cases, neutrino oscillation can be used as an independent check in the event that 
such new physics is discovered at the LHC 

All the other models are already well constrained by existing experiments, either indirectly 
by low-energy precision measurements, or by direct searches at colliders. Generically, the 
couplings and masses of new particles that couple only to leptons are well constrained by 
lepton universality, while their contribution to neutrino oscillation tend to be suppressed 
since they only interact with the electrons in matter. This tends to render the existing 
bound stronger than the potential bound from |^| > ^o- 

Topcolor assisted technicolor, and R-parity violating LQD couplings involve interactions 
with the quarks in matter, but they too belong to the list of already well-constrained models. 
For the Z' in topcolor assisted technicolor, the proton and electron contributions to neutrino 
oscillation cancel, just as for the Standard Model Z, and the coupling is also fixed to a small 
value, which results in a weak bound from |,^| < ^o- For the LQD coupling, restriction to 
minimal supergravity provided an extra constraint which strengthened the existing bound. 

The fact that only a limited number of models (at least among those we surveyed) can be 
well constrained by |^| < means, conversely, that if a non-zero ^ is observed in neutrino 
oscillation, the list of possible new physics that could lead to such an effect is also hmited. 
This could, in principle, help distinguish among possible new physics which have the same 
type of signature {e.g. a leptoquark which may, or may not be generation diagonal) at the 
LHC. 
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